PROOF OF THE MONOTONE COLUMN 
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In memoriam Julius Borcea 

Abstract. Let A be an n-by-n matrix of real numbers which are weakly de- 
creasing down each column, Z„ = diag(zi, . . . , z n ) a diagonal matrix of inde- 
tcrminates, and J n the rt-by-n matrix of all ones. We prove that per( J n Z n +A) 
is stable in the Zj, resolving a recent conjecture of Haglund and Visontai. This 
immediately implies that per(2j n + A) is a polynomial in z with only real 
roots, an open conjecture of Haglund, Ono, and Wagner from 1999. Other 
applications include a multivariate stable Eulerian polynomial, a new proof of 
Grace's apolarity theorem and new permanental inequalities. 



1. The monotone column permanent conjecture. 

Recall that the permanent of an n-by-n matrix H = (/%) is the unsigned variant 
of its determinant: 

n 

per(F) = ^2 Yi h iMi)> 

iE6„ t=l 

with the sum over all permutations a in the symmetric group & n . A monotone 
column matrix A = (ay ) has real entries which are weakly decreasing reading down 
each column: that is, ay > a%+i,j for all 1 < i < rt — 1 and 1 < j ' < n. Let J„ be 
the n-by-n matrix in which every entry is 1. 

The Monotone Column Permanent Conjecture (MCPC). 

If A is an n-by-n monotone column matrix then per(z J n + A) is a polynomial in z 
which has only real roots. 

The MCPC first appears as Conjecture 2 of [7]. (Originally with increasing 
columns - but the convention of decreasing columns is clearly equivalent and will be 
more natural later.) Theorem 1 of [7] proves the MCPC for monotone column ma- 
trices in which every entry is either or 1. Other special cases appear in [SJ [71 [51 [H] , 
either for n < 4 or for rather restrictive conditions on the entries of A. In this 
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paper we prove the MCPC in full generalitjo In fact we prove more. The following 
multivariate version of the MCPC originates in [8]. Let Z n = diag(zi, . . . , z n ) be 
an n-by-n diagonal matrix of algebraically independent commuting indeterminates 
z = {zi, . . . , z n J. A polynomial f(zi, . . . , z n ) £ C[z] is stable provided that either 
/ = or whenever Wj £ C are such that lm(wj) > for all 1 < j < n, then 
f(wi, . . . , w n ) ^ 0. A stable polynomial with real coefficients is real stable. 

The Multivariate MCP Conjecture (MMCPC). 

If A is an n-by-n monotone column matrix then per( J n Z n + A) is a real stable 
polynomial in M[z]. 

Note that J n Z n + A = (zj + dij), so that Zj is associated with the j-th column, 
for each 1 < j < n. We also write per( J n Z n + A) = pcr(zy +aij) as it seems clearer. 
The MMCPC implies the MCPC since if one sets all Zj = z, then per(z + a^) is 
a polynomial in one variable with real coefficients; this diagonalization preserves 
stability (see Lemma 12.1( c) below), and a univariate real polynomial is stable if 
and only if it has only real roots. 

In Section 2 we review some results from the theory of stable polynomials which 
are required for our proofs. In Section 3 we reduce the MMCPC to the case of {0, 1}- 
matrices which are weakly decreasing down columns and weakly increasing from 
left to right across rows - these we call Ferrers matrices for convenience. Then we 
further transform the MMCPC for Ferrers matrices, derive a differential recurrence 
relation for the resulting polynomials, and use this and the results of Section 2 to 
prove the conjecture by induction. In Section 4 we extend these results to sub- 
permanents of rectangular matrices, derive a cycle-counting extension of one of 
them, discuss a multivariate stable generalization of Eulerian polynomials, present 
a new proof of Grace's apolarity theorem and derive new permanental inequalities. 

2. Stable polynomials. 

Over a series of several papers, Borcea and Branden have developed the theory 
of stable polynomials into a powerful and flexible technique. The results we need 
are taken from [H [21 H] ; see also Sections 2 and 3 of [T3] ■ 

Let "K — {w € C : Im(u>) > 0} denote the open upper half of the complex 
plane, and let 3~C denote its closure in C. As above, z = {z%, . . . , z n } is a set of n 
indeterminates. For / £ C[z] and 1 < j < n, let deg z .(/) denote the degree of Zj 
in /. 

Lemma 2.1 (see Lemma 2.4 of [13] ). These operations preserve stability of poly- 
nomials in C[z] . 

(a) Permutation: for any permutation a £ 6,„ / M> f(z a i^)^ ■ ■ ■ , z a ^). 

(b) Scaling: for c £ C and a £ 1™ with a > 0, / H> c/(aizi, . . . , a n z n ). 

(c) Diagonalization: for 1 < i < j < n, f i-> f(z)\ Zi=Zj . 

(d) Specialization: for a £ Di, f M> f(a, z%, . . . , z n ). 

(c) Inversion: if deg zi (f) = d, f i-> zff(-z^ 1 ,z 2 , . . .,z n ). 

(f) Translation: / \-> fi = f{z\ + t, z 2l ■ ■ ■ , z n ) £ C[z, t}. 

(g) Differentiation: / n> df(z)/dzi. 



Thus, it was not really a permanent conjecture. 
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Proof. Only part (f ) is not made explicit in [TJ [21 [13] . But clearly if z\ G Jt and 
i £ M then z\ + t G Jf , from which the result follows. □ 

Of course, parts (d,e,f,g) apply to any index j as well, by permutation. Part (g) 
is the only difficult one - it is essentially the Gauss-Lucas Theorem. 

Lemma 2.2. Let f(z,t) e C[z, i] be stable, and let 

d 

with /d(z) ^ 0. Then /fe(z) is stable for all < k < d — deg t (/). 

Proof. Consider anyO < k < d. Clearly /^(z) is a constant multiple of d k f(z, t)/dt k \ t= o, 
which is stable by Lemma 12. II (c. g). □ 

Polynomials g, h £ K[z] are in proper position, denoted by g <C ft, if the poly- 
nomial ft + ig is stable. This is the multivariate analogue of interlacing roots for 
univariate polynomials with only real roots. 

Proposition 2.3 (Lemma 2.8 of [2 and Theorem 1.6 of pQ). Let g,h€ K[z]. 
(a) Then h <C g if and only if g + th G R[z, t] is stable. 

(a) Then ag + bh is stable for all a, b G K i/ and on/?/ i/ either h <C g or g ^ h. 

It then follows from Lemma I2.1f d,g) that if ft. 5 then both ft, and 17 are stable 
(or identically zero). 

Proposition 2.4 (Lemma 2.6 of pQ). Suppose that g £ K[z] is stable. Then the 
sets 

{ft S M[z] : g < ft} and {ft £ E[z] : ft < 5} 
are convex cones containing g. 

Proposition 2.5. Lei fee a reaZ vector space, <fi : V n R a multilinear form, 
and ei, . . . , e„, V2, ■ ■ ■ , v n fixed vectors in V . Suppose that the polynomial 

4>(ei,v 2 + z 2 e 2l ...,v n + z n e n ) 

in M[z] is not identically zero. Then the set of all V\ G V for which the polynomial 

<j>(vi + zid,v 2 + z 2 e 2 , ...,v n + z n e n ) 

is stable is either empty or a convex cone (with apex 0) containing e± and —e\. 

Proof. Let C be the set of all v\ 6 V for which the polynomial <p(vi + z\e\, v 2 + 
z 2 e 2 , ...,«„ + z n e n ) is stable. For v € V let F v = <j>(v, v 2 + z 2 e 2 , ...,«„ + z n e n ). 
Since 

4>{vi + ziei,v 2 + z 2 e 2 , ...,v n + z n e n ) = F Vl + z x F ei , 

we have C = {v E V : F ei « Moreover since F\ v+f2W = XF V + ^,F W it follows 
from Proposition 12.41 that C is a convex cone provided that C is non-empty. If C 
is nonempty then F v + Z\F ex is stable for some v £ V. But then F ei is stable, and 
so is 

(±1 + Zi)F ei = cj>{±ei + ziei,v 2 + z 2 e 2 , ...,v n + z n e n ) 
which proves that ±ei G C. □ 
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(Of course, by permuting the indices Proposition 12.51 applies to any index j as 
well.) 

Let C[z] ma denote the vector subspace of multiaffine polynomials: that is, poly- 
nomials of degree at most one in each indeterminate. 

Proposition 2.6 (Theorem 5.6 of j4]). Let f G R[z] ma be multiaffine. The following 
are equivalent: 

(a) / is real stable. 

(b) For all 1 < i < j < n and all a G E™, / l (a)/ J (a) — /(a)/y(a) > 0, in which the 
subscripts denote partial differentiation. 

A linear transformation T : C[z] — > C[z] preserves stability if T(/)(z) is stable 
whenever /(z) is stable. 

Proposition 2.7 (Theorem 1.1 of [2]). Let T : C[z] ma -)• C[z] be a linear transfor- 
mation. Then T preserves stability if and only if either 

(a) T(f) = rj{f) ■ p for some linear functional 7/ : C[z] ma — > C and stable p G C[z], 
or 

(b) i/ie polynomial Gt(z, w) = Trj™ =1 (z.j + Wj) *s stable in C[z,w]. 

The complexification Tfc : C[z] — > C[z] of a linear transformation T : R[z] — > R[z] 
is dehned as follows. For any / G C[z] write / = g + ih with g,h £ R[z], which can 
be done uniquely. Then T c (f) = T{g) + iT(h). Let also Tfc = T(g) - iT(/i). 

Proposition 2.8. Let T : M[z] — > M[z] be a linear transformation and let f G C\z] 
be a stable polynomial. If T preserves real stability then either 7c(/) is stable, or 
Tc{.f) ^ stable. 

Proof. Let / = g + ih G C[z] with g,h £ K[z], and assume that / is stable. Then 
h -c 5 by definition. By Proposition 12 ,3f b) . it follows that a/i + bg is real stable for 
all a, b G M. Therefore, aT(h) + bT(g) is real stable for all a, 6 G R. By Proposition 
HJ2» again, either T(/i) < T( 5 ) or T( 5 ) < T(h). By Proposition E^a), either 
T(g) + iT(/i) = T c (/) is stable, or T(g) - \T(h) = f c (f) is stable. □ 

3. Proof of the MMCPC. 
3.1. Reduction to Ferrers matrices. 

Lemma 3.1. If per(zj + ay) is stable for all Ferrers matrices, then the MMCPC 
is true. 

Proof. If per(zj + a^) is stable for all Ferrers matrices, then by permuting the 
columns of such a matrix, the same is true for all monotone column {0, l}-matrices. 
Now let A = iflij) be an arbitrary n-by-n monotone column matrix. We will show 
that per(2j + ay) is stable by n applications of Proposition ^. 51 

Let V be the vector space of column vectors of length n. The multilinear form 
4> we consider is the permanent of an n-by-n matrix obtained by concatenating n 
vectors in V. Let each of e\, . . . , e„ be the all-ones vector in V. 

Initially, let vi,V2, ■ ■ ■ , v n be arbitrary monotone {0, 1}- vectors. Then 4>{v\ + 
z\e\, ...,v n + z n e n ) = per(J n Z n + H) for some monotone column {0, l}-matrix H. 
One can specialize any number of Vj to the zero vector, and any number of Zj to 
1, and the result is not identically zero. By hypothesis, all these polynomials are 
stable. 
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Now we proceed by induction. Assume that if V\, . . . ,Vj-i are the first j — 
1 columns of A, and if Vj,...,v n are arbitrary monotone {0, l}-columns, then 
4>(vi + ziei, ...,v n + z n e n ) is stable. (The base case, j — 1, is the previous para- 
graph.) Putting Vj — and Zj = 1, the resulting polynomial is not identically 
zero. By Proposition 12.51 (applied to index j), the set of vectors Vj such that 
4>(vi + z±ei, . . . , v n + z n e n ) is stable is a convex cone containing iej. Moreover, it 
contains all monotone {0, l}-columns, by hypothesis. Now, any monotone column of 
real numbers can be written as a nonnegative linear combination of —e\ and mono- 
tone {0, l}-columns, and hence is in this cone. Thus, we may take v\ , . . . , , Vj to 
be the first j columns of A, Wj+i, . . . ,v n to be arbitrary monotone {0, l}-columns, 
and the resulting polynomial is stable. This completes the induction step. 

After the n-th step we find that per( J n Z n + A) is stable. □ 

3.2. A more symmetrical problem. Let A = (ay) be an n-by-n Ferrers matrix, 
and let z = {zi, . . . ,z n }. For each 1 < j < n, let y.j = (zj + l)/zj, and let 
Y n = diag(yi, . . . , y n ). The matrix obtained from J n Z n + A by factoring Zj out of 
column j for all 1 < j < n is AY n + J n — A = (aijyj + 1 — oy). It follows that 

per(zj + a,ij) = z x z 2 ■ ■ ■ z n ■ per(ayj/j + 1 - fly). (3.1) 

Lemma 3.2. For a Ferrers matrix A — (a^-), per(zj + a^) is stable if and only if 
per iflijyj + 1 — ay) is stable. 

Proof. The polynomials are not identically zero. Notice that Im(zj) > if and only 
if Im(yj) = Im(l + zj 1 ) < 0. If per(z J + ay) is stable, then per(ayj/^ + 1 — ay) ^ 
whenever lm(yj) < for all 1 < j <n. Since this polynomial has real coefficients, 
it follows that it is stable. The converse is similar. □ 



The set of n-by-n Ferrers matrices has the following duality A h- > A y = J n — A T : 
transpose and exchange zeros and ones. That is, A v = (dy) in which ay = 1 — ciji 
for all 1 < i,j < n. Note that (^4 V ) V = A. However, the form of the expression 
per(ayyj + 1 — ay) is not preserved by this symmetry. To remedy this defect, 
introduce new indeterminates x = {xi, . . . ,x n } and consider the matrix B(A) = 



(bij) with entries bi 



jyj + (1 — aij)xi for all I < i, j < n. For example, if 



.4 



1111 

111 

111 

1 1 





then B(A) 



xi y 2 V3 V5 

X2 x 2 y 3 2/4 V5 

%3 x 3 y 3 j/4 y 5 

X 4 X 4 X 4 2/4 2/5 

^5 ^5 ^5 



For emphasis, we may write B(A;x;y) to indicate that the row variables are 
x and the column variables are y. The matrices B(A) and B(A V ) have the same 
general form, and in fact 



per(B(A v ;x;y)) = pev(B(A; y; x)). 



(3.2) 



Clearly peT(B(A)) specializes to per(ayt/j + 1 — ay) by setting Xi = 1 for all 



1 < i < n. We will show that per(_B(/l)) is stable, for any Ferrers matrix A. 
Lemmas 0;d),[321 andEU this will imply the MMCPC. 



By 
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□ 



□ 



□ 



□ 



□ 



□ 



□ 



□ 



□ 



□ 



□ 



Figure 1. ct=316524 maps to ir(a) = 3 1 4 5 2. 



3.3. A differential recurrence relation. Next, we derive a differential recur- 
rence relation for polynomials of the form per(B(A)), for A an n-by-n Ferrers 
matrix. There are two cases: either a nn — or a nn = 1. Replacing A by A^ and 
using (|3.2[) . if necessary, we can assume that a nn = 0. 

Lemma 3.3. Let A = (oy) be an n-by-n Ferrers matrix with a nn — 0, let k > 1 be 
the number of 's in the last column of A, and let A be the matrix obtained from 
A by deleting the last column and the last row of A. Then 

per(B(A)) = kx n per(B(A°)) + x n y n dpci(B(A°)), 

in which 



71 — k 71 — 1 

d = V — + V — 



Proof. In the permutation expansion of per(B(A)) there are two types of terms: 
those that do not contain y n and those that do. Let T a be the term of per(B(A)) 
indexed by a £ & n . For each n — k + 1 < i < n, let Gi be the set of those terms To- 
such that o~(i) = n; for a term in 6^ the variable chosen in the last column is Xi. 
Let D be the set of all other terms; for a term in D the variable chosen in the last 
column is y n . 

For every permutation a G ©„, let (i a ,ja) be such that o~(i a ) — n and cr(n) = j a , 
and define ir(a) G & n -i by putting 7r(i) = o~(i) ifi=£ i a , and 7r(i CT ) = j a (if i a ^ n). 
Let Tn( a ) be the corresponding term of pei(B(A°)). See Figure 1 for an example. 
Informally, 7r(cr) is obtained from <r, in word notation, by replacing the largest 
element with the last element, unless the largest element is last, in which case it is 
deleted. 

For each n — k + 1 < i < n, consider all permutations a indexing terms in 6^. 
The mapping T a M> T n / a s is a bijection from the terms in Gi to all the terms in 
per(B(A°)). Also, for each a e C^, T a — x n T^^ a y Thus, for each n — k + 1 < i < n, 
the sum of all terms in Gi is x n per(B(A )). 

Next, consider all permutations a indexing terms in D. The mapping T CT i— > 
T-nta) is (n — fc)-to-one from D to the set of all terms in per(B(A°)), since one 
needs both ir{a) and i a to recover a. Let v a be the variable in position {i a ,j a ) 
of B(A°). Then v a T a = XnynT^^y It follows that for any variable w in the set 
{xi, . . . , x n -k, j/i, ... , y n -i}, the sum over all terms in D such that v a — w is 

d 



x n y n —per(B(A )). 



Since v a is any element of the set {xi, . . . , x r , 
sum of all terms in D is x n y n dper(B(A°)). 



-k,Vi, 



,y n -i}, it follows that the 
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The preceding paragraphs imply the stated formula. □ 
3.4. Finally, proof of the MMCPC. 

Theorem 3.4. For any n-by-n Ferrers matrix A, ~per(B(A)) is stable. 

Proof. As above, by replacing A by A v if necessary, we may assume that a\ n = 
0. We proceed by induction on n, the base case n = 1 being trivial. For the 
induction step, let A be as in Lemma 13.31 By induction, we may assume that 
pei(B(A°)) is stable; clearly this polynomial is multiaffine. Thus, by Lemma [3~3l it 
suffices to prove that the linear transformation T = k + y n d maps stable multiaffine 
polynomials to stable polynomials if k > 1. This operator has the form T — 
k + z m 5Zj=i d/dzj (renaming the variables suitably). By Proposition ^ . 71 it suffices 
to check that the polynomial 

G T (z,w) = T JJ(2j +Wj) 



3=1 3 3 J 3=1 

is stable. If zj and Wj have positive imaginary parts for all 1 < j < m then 

1 m~\ 

e = - + E — 

Zm i 
3 = 1 



has negative imaginary part (since k > 0). Thus z m £ ^ 0. Also, Zj + Wj has 
positive imaginary part, so that Zj + Wj ^ for each 1 < j < m. It follows 
that Gt(z,w) 7^ 0, so that Gt is stable, completing the induction step and the 
proof. □ 

Proof of the MMCPC. Let A be any n-by-n Ferrers matrix. By Theorem 13.41 
pei(B(A)) is stable. Specializing Xi = 1 for all 1 < i < n, Lemma 12.1( d) im- 
plies that per(ay2/j + 1 — a%j) is stable. Now Lemma 13.21 implies that per(zj + ay) 
is stable. Finally, Lemma \3 . 1 1 implies that the MMCPC is true. □ 



Henceforth, we shall refer to the MMCPT - "T" for "Theorem" 



4. Further results. 

4.1. Generalization to rectangular matrices. We can generalize Theorem 13.41 
to rectangular matrices, as follows. Let A = (a^ ) be an m-by-n matrix. As in the 
square case, A is a Ferrers matrix if it is a {0, l}-matrix that is weakly decreasing 
down columns and increasing across rows. The matrix B(A) is constructed just 
as in the square case: B(A) = (bij) = {a^yj + (1 — aij)xi), using row variables 
x = {x\, . . . , x m } and column variables y = {j/i, . . . , y n }. For emphasis, we may 
write B(A; x; y) to indicate that the row variables are x and the column variables 
arc y. The symmetry A i— > A y takes an m-by-n Ferrers matrix to an n-by-m Ferrers 
matrix. 

Now let k < min{m, n}. The k-permanent of an m-by-n matrix H — (hij) is 

R C i£R 



s 
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in which R ranges over all fc-element subsets of {1, . . . , m}, C ranges over all k- 
element subsets of {1, . . . , n}, and (3 ranges over all bijections from R to C. (Note 
that per (i?) = 1 for any matrix H.) In the case k = m = n this reduces to the 
permanent of a square matrix. 

For an m-by-n Ferrers matrix A and k < min{m, n), note that 

per fc (i?(A v ;x; y )) = per fc (B(A; y; x)). (4.1) 

Thus, replacing A by A if necessary, we may assume that m < n. 

Proposition 4.1. Let A — (ay) &e an m-by-n Ferrers matrix. Then per k (B(A)) 
is stable for all k < min{m, n}. 

Proof. Using (14. ip . if necessary, we may assume that m < n. 

We begin by showing that per m (B(A)) is stable. Let A' be the ro-by-n Ferrers 
matrix obtained by concatenating n~m rows of O's to the bottom of A. One checks 
that 

per(B(A')) = (n - to)! x n x n -i ■ ■ ■ x m+1 ■ pex m (B(A)). 
By Theorem 13. A\ per(B(A')) is stable, and it follows easily that per m (B(A)) is 
stable. 

Now, let J m ,n be the m-by-n matrix of all l's. Then 

pex m (B(A) + U m . n ) = per k (B(A)) ( U _ J (m - k)lt m ' k . 

By Lemma [2.ir c.f). this polynomial is stable. Extracting the coefficient of f n ~ k 
from this, and dividing by ( n _ k )( m — Lemma [2.21 shows that per fc (B(^4)) is 
stable for all < k < m. □ 

Proposition [4j] suggests the idea of a similar generalization of the MMCPT: is it 
true that per fc ( J m , n Z n + A) is stable for every m-by-n monotone column matrix A 
and k < min{m, n}? This conjecture originates in [8]. One cannot derive this from 
Proposition 14.11 however, because there is no analogue of (|3.1j) for fc-permanents. 
Nonetheless, we can prove this result for half the cases. 

Proposition 4.2. Let A be an m-by-n monotone column matrix with to > n, and 
let k < n. Then per k (J m ^ n Z n + A) is stable. 

Proof. Let A' be the m-by-m matrix obtained from A by concatenating to — n zero 
columns to the right of A. Then 

per(J m Z m + A') = (to - n)\z m z m -i ■ ■ ■ z n+ i ■ per n ( J m , n Z n + A). 

Since pei(J m Z m + A') is stable, it follows that pe£ n (J m , n Z n + A) is stable. By 
Lemma [2~l7 c.f). it follows that per n (J m . n Z n + A + tJ m , n ) is stable. Extracting the 
coefficient of t n ~ k from this, and dividing by ( m Zfe)(^ — Lemma [2~2l shows that 
peT k (J m ,nZ n + A) is stable for all < k < n. □ 

4.2. A cycle-counting generalization. Theorem 13.41 can be generalized in an- 
other direction, as follows. 

For each permutation a € & n , let cyc(<r) denote the number of cycles of a. For 
an indeterminate a and an n-by-n matrix H = (/i,j), the a-permanent of H is 

n 

per(if;a)= £ a^JJ^-W- 
CT ee„ i=i 
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The numbers cyc(er) behave well with respect to the duality A H> A v : for any 
Ferrers matrix A, 

per(£(A v ; x; y); a) = per(B(A; y; x); a). (4.2) 

Lemma 4.3. Let A = (oj 3 -) be an n-by-n Ferrers matrix with a nn = 0, let k > 1 be 
the number of 's in the last column of A, and let A° be the matrix obtained from 
A by deleting the last column and the last row of A. Then 

per(B(A); a) = (a + k — l)x n per(B(A°); a) + x n y n dpei(B(A°); a), 

with d as in Lemma\3.3\ 



Proof. Adopt the notation of the proof of Lemma 13731 To obtain the present result, 
observe that if T a is in S n then cyc(cr) = 1 + cyc(7r(cr)), and otherwise cyc(cr) = 
cyc(7r(V)). □ 

Proposition 4.4. For a > and A a Ferrers matrix, pev(B(A); a) is stable. 

Proof. Reprising the proof of Theorem 13.41 it suffices to show that an operator of 
the form 

m—l 

T=(a + k-l) + z n Y^ d / dz * 

i=l 

preserves stability when k > 1 . The argument of the proof of Theorem 13.41 works 
when a > 0. □ 

For a > and A a Ferrers matrix, specialize all X{ = 1 and diagonalize all yj = z 
in per(£?(j4); a). By Lemma I2.1f c.d). the result is a (univariate) polynomial with 
only real roots. This special case is also implied by Theorem 2.5 of [6]. 

4.3. Multivariate stable Eulerian polynomials. Given a permutation a £ ©„, 
viewed as a linear sequence <t(1)<j(2) • • ■ <r(n), let L(er) denote the result of the 
following procedure. First form the two-line array 

1 2 ••• n 
v cr(l) <j(2) ••• a(n) / 

Then, viewing T(er) as a map sending i to o~(i), break T(a) into cycles, with the 
smallest element of each cycle at the end, and cycles listed left-to-right with smallest 
elements increasing. Finally, erase the parentheses delimiting the cycles to form a 
new linear sequence L(<r). For example, if a = 341526978, then T(cr) has cycle 
decomposition (31) (452) (6) (987) and L(a) = 314526987. 

Let P(<t) denote the placement of n nonattacking rooks on the squares (i,o~(i)), 
1 < i < n. As noted by Riordan 10 , rooks in P(<r) that occur above the diagonal 
correspond to exceedences (values of i for which o~(i) > i) in a and descents (values 
of i for which a{i) > a(i + 1)) in L(cr). Hence 

J2 g"**.') = z dos(CT) , (4.3) 

where exc(cr) is the number of exceedences and des(cr) the number of descents of a. 

The polynomial in (|4.3I) is known as the Eulerian polynomial. It is one of the 
classic examples in combinatorics of a polynomial with only real roots. Let E n — 
(eij) denote the n-by-n matrix with = if i > j and = 1 if i < j. Then 
pei(B(E n ; x; y)) is stable by Theorem l3.4l Let 1 be the vector (of appropriate size) 
of all ones. From the above discussion, per(i?(i?„; 1; z, . . . , z)) equals the Eulerian 



T(a) = 
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polynomial; by Lemma I2.1f c.d) this is a univariate stable polynomial with real 
coefficients, so it has only real roots. Similarly, we see that 

per(B(£„;l;y))= £ J[ y a(l) (4.4) 

<ree„ <r(i)>tj(i+l) 

= E II ( 4 - 5 ) 

crE@„ a(i)>i 

is stable in {yi, 1/2, • ■ • , 2/n} (but ?/i does not really occur). 

Letting / = pei(B(E n ; l;y)), note that the partial of / with respect to j/j, 
evaluated at all y- variables equal to 1 , equals the number of permutations in & n 
that have i as a "descent top" , i.e. have the property that i is followed immediately 
by something less than i. Denoting this number by Top(i; n), applying Proposition 
12.61 to / we get 

Top(i;n)Top(j;n) > n\Top(i,j;n), (4.6) 

where Top(i, j;n) is the number of permutations in 6 n having both i and j as 
descent tops, with 2 < i < j < n. Dividing both sides of the above equation by 
n! 2 shows that occurrences of descent tops in a uniformly random permutation are 
negatively correlated. 

More general forms of (|4.4|) can be defined which still maintain stability. First 
of all, cycles in T(cr) clearly translate into left-to-right minima in L(cr), and so by 
Proposition 14.41 the polynomial 

e « LRmin(CT) n few 

itE6„ cr(i)>a(i+l) 

is stable in the y for a > 0, with LRmin(cr) denoting the number of left-to-right 
minima of a. 

Secondly, the sum in (|4.4|) can be replaced by a sum over permutations of a 
multiset. For a given vector v = (vx,...,v t ) & N*, let N(v) = {l Vl 2 V2 ■ ■ ■ t Vt } 
denote the multiset with Uj copies of i, and M(v) the set of multiset permutations 
of N(v) (so for example M(2, 1) = {112,121,211}). Riordan [TD] noted that if 
we map our previous sequence L(a) to a multiset permutation m(cr) by replacing 
numbers 1 through v\ in L(a) by l's, numbers v\ + 1 through Ui + «2 by 2's, etc., 
we get a Y[ w<! to 1 map, and furthermore certain squares above the diagonal where 
rooks in P(cr) correspond to descents in L{o) no longer correspond to descents in 
m(cr). For example, if v% — 2, then 1 and 2 in L(a) both get mapped to 1 in m(a), 
so a rook on square (1, 2) no longer corresponds to a descent. 

Let n denote the sum of the coordinates of v, and let Y(v) be the sequence of 
variables obtained by starting with yi, . . . , y n and setting the first v\ y- variables 
equal to yi, the next w 2 y- variables equal to 2/2, etc. Then if E(v) is the Ferrers 
matrix whose first v\ columns are all zeros, the next V2 columns have ones in the 
top v\ rows and zeros below, the next V3 columns have ones in the top v\ + vi rows 
and zeros below, etc., an easy extension of the argument above implies 

(l/n^)per(5(£(v;l;Y(v))))= ]T JJ f'W ( 4J ) 

i crGM(v) <r(i)>a{i+l) 

is stable in the yi. This contains Simion's result [11] . that the multiset Eulerian 
polynomial has only real roots. If v has all ones, i.e. iV(v) is a set, it reduces to our 
previous result. Finally, we note that this argument also shows that if we replace 
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the condition o~(i) > o~(i + 1) by the more general condition a(i) > o~(i + 1) +j — 1, 
for any fixed positive integer j, we still get stability. 

4.4. Grace's Apolarity Theorem. Univariate complex polynomials f(t) = X)fe=o 
and g(t) = X)fe=o Ck)^ k are a P°^ ar if a nb n ^ and 



n 

^ \k 



{-l) n - k a k b n - k = 0. 



Lemma 4.5. Let /(f) = X)fe=o (fc) afe ^ fc arl ^ = Sfe=o (fc)^ fc ^ e com pl ex poly- 
nomials of degree n. Let the roots of f(t) be z\, . . . , z n and let the roots of g(t) be 
u>i,...,w n - Then 



n 

\ - / n 



^ \k 

k=0 



(-1)" k a k b n - k = n\a n b n peT(w l - Zj). 



In particular, f(t) and g(t) are apolar if and only ifper(wi — Zj) = 0. 

Proof. It suffices to prove this for monic polynomials f(t) and g(t). For each per- 
mutation a € & n there are 2™ terms in per(wi — Zj), since for each 1 < i < n cither 
Wi or — z CT (j) can be chosen. For each subset R of rows of size k, and subset C of 
columns of size n — k, the monomial w R z is produced kl(n — k)l times. Since 
(— l) k a n - k is the fc-th elementary symmetric function of {z\, . . . , z n }, and similarly 
for (— l) k b n _ k and {tui, . . . ,w n }, the result follows. □ 

Lemma 4.6. Let f(t) and g(t) be polynomials of degree n. Let t t— > <j)(t) — (at + 
b)/(ct + d) be a Mobius transformation, with inverse — (at + j3)/( r )t + S). 

Let f(t) = (-ft + tf)"/(0 -1 (t)) and g(t) = (jt + 5) n g(<j)- x (t)) have degree n. Then 
f(t) and g(t) are apolar if and only if f(t) and g(t) are apolar. 

Proof. Let the roots of f(t) be z\, . . . , z n and let the roots of g(t) be w\, . . . , w n . 
Then the roots of f(t) are <fi(zi), . . . , <f>(z n ) and the roots of g(t) are (f>(wi), . . . , (j)(wn)- 
Consider the permanent per(0(i«i) — <p(zj)). The (i, j)-entry of this matrix is 

awi + b azj + b (ad — bc)(wi — Zj) 
cwi + d czj + d (cu>i + d) (czj + d) 

Factor (cwi + c?) _1 out of row i, factor (czj + d)^ 1 out of column j, and factor 
ad — be out of every row. Therefore 

pei(<j>(wi) - (j>(zj)) = ^ — ^ h ^ —- ■ per(wj - zj). 

Uh=i ( cw h + d)(cz h + d) 

Since the prefactor on the right-hand side is neither zero nor infinite, the result 
follows from Lemma WM Q 



A circular region A is a proper subset of C that is either open or closed, and is 
bounded by either a circle or a straight line. 

Theorem 4.7 (Grace's Apolarity Theorem). Let f(t) and g(t) be apolar polyno- 
mials. If every root of g(t) is in a circular region A, then f(t) has at least one root 
in A. 
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Proof. It suffices to prove this for monic polynomials f(t) and g(t), and for open 
circular regions A since a closed circular region is the intersection of all open circular 
regions which contain it. Let t i— > <fi(t) be a Mobius transformation that maps A to 
the upper half-plane !K. By Lemma 14.61 it suffices to prove this when the circular 
region is 5f itself. Let the roots of f(t) be z\, . . . , z„ and let the roots of g(t) be 
wi, ...,w„. 

If all of Wi, . . . , w n are real numbers, then by permuting the rows of + Zj) 
we can assume that w\ < • • • < w n without changing the value of pei(wi + Zj). By 
the MMCPT, per(u>i + Zj) is a real stable polynomial in z%, . . . , z n . In other words, 
the transformation T : R[z] — »• M[z] dchncd by 

T (f( z )) = P<2i(w l +Zj), 

where /(z) = (z + w±)(z + W2) • ■ ■ (z + w„) preserves real stability. This is a 
linear transformation. Suppose that f(z) G C[z] is stable. By Proposition 12.81 
either Tc(/(z)) is stable or Tc(f(z)) is stable. Diagonalizing by setting Zj = z for 
ah 1 < j < we see that Tc(f)(z, . . . , z) = nlf(z), so that Tc(/(z)) is stable. 
Therefore per(i«i + Zj) is a stable polynomial in C[z], for all wi, . . . ,w n G 5£. 
Therefore per(i«i + Zj) is a stable polynomial in C[w,z]. Actually it satisfies a 
stronger stability property. Namely if Zj G 'K and G IK for all 1 < j < n, then 
per(wi + Zj) ^ 0. Indeed if we fix Q G "X for all 1 < j < n, then the polynomial 
per(wi + Cj) G C[w] is cither identically zero or stable (and not identically zero) by 
Lemma |2. II (d). Clearly per(wi + Q) is not identically zero. 

Now assume that Wi G IK for all 1 < i < n. Arguing for a contradiction, 
assume that Zj IK for all 1 < j < n. Then — Zj G K for all 1 < j < n. 
Hence per(wi — Zj) ^= by the argument given in the previous paragraph. But 
per(u>i — Zj) — by Lemma 14.51 since f(t) and g(t) are apolar. This contradiction 
completes the proof. □ 

Our proof of Grace's apolarity theorem relies on MMCPT. To avoid going in 
circles we should ensure ourselves that the proof of MMCPT does not use any form 
of Grace's apolarity theorem. In fact what we use in the proof of MMCPT is that 
condition (b) in Proposition 12.71 implies that the operator preserves stability, the 
proof of which does not use Grace's apolarity theorem. 

4.5. Permanental inequalities. If A is a n-by-n matrix and SC[n] = {l,...,n} 
let As be the matrix obtained by replacing the columns indexed by S by columns 
of all ones. 

Corollary 4.8. Let A be an n-by-n monotone matrix, S C [n], and i,j G [n] \ S. 

Then 

per(A Su{i} )per(A 5u{j} ) > per(A Su{i j } )pei(As). 

Proof. The proof follows immediately from Proposition 12.61 applied to per(zj + 
fly). □ 

Note that Corollary 14.81 can be interpreted as that the permanent is pairwise 
negatively associated in columns. 

The generating polynomial of a discrete measure /1 : 2^ —¥ K>o is given by 

g„( Z )= ^)Y[ Zj . 

S62W j6S 
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The measure /i is Rayleigh if 

^(x)^(x) > |!^L(x)G AI (x), (4.8) 

ozi ozj oziOz-j 

for all x G R> , and it is called strongly Rayleigh if (|4.8j) holds for all xel". We 
refer to [3l [12] for more information on Rayleigh and strongly Rayleigh measures. 

Suppose that A is an n-by-n monotone matrix with nonnegative coefficients. 
Consider the discrete measure fiA ■ 2^ M>o defined by ha(S) = per(^s). 
By Proposition 12.61 we see that ha is strongly Rayleigh. This fact entails many 
inequalities. 

Corollary 4.9. Suppose that A is an n-by-n monotone matrix with nonnegative 
coefficients. Then 



for all S,TC [n] . 
Moreover 



per(A s )per(^ T ) > per(^ SuT )per(ylsnT), 



per(A) < si ••■ s n — - 



where Si is the sum of the elements in the ith column. 

Proof. The first inequality holds for all Rayleigh measures, see [12j Theorem 4.4]. 

Let n(S) = per(A[„]\ s )/n!. By the above h(S)h(T) > fi(S U T)h(S n T) for all 
S,T C [n]. Thus n(S U T) < fj,(S)K T ) whenever S n T = 0, and after iteration 
^(M) < M({1}) ' • ' Mi 71 -})- The proof now follows by observing that = Si/n 

for all i € [n] . 

One can also prove the last inequality by an elementary argument. If there are 
two different consecutive elements b > a in a column of A, replace these by their 
average to obtain the matrix A' . It is plain to see that per(A) < per(A'). Iterating 
this procedure it follows that per (A) < per(i3), where each element in column i of 
B is equal to Si/n. The inequality now follows. □ 



The second inequality in Corollary 14.91 can be compared with the Van der Waer- 
den conjecture which asserts that the permanent of a doubly stochastic matrix is 
greater than or equal to n\/n n . The Van der Waerden conjecture which was stated 
in 1926 was proved by Falikman in 1981; the case of equality was proved by Ego- 
rychev. Gurvits has recently provided a beautiful proof of a vast generalization of 
the Van der Waerden conjecture using stable polynomials, see (5] and [T3l Section 
8]- 

If / : 2 1"! -> K, and ji is a discrete measure on 2'™' we let 



f fdfJ- = J2 /(SMS)- 



A measure fj- on 2'™' is negatively associated if for all increasing functions /, g 
2M — >. R depending on disjoint sets of variables 



fgd^i dfi< fdfi / gdfi 



see e.g. [3]. 
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Corollary 4.10. Suppose that A is an n-by-n monotone matrix with nonnegative 
coefficients. Then the discrete measure \ia ■ 2W -> R> , defined by /j, A (S) = 
pev(As), is negatively associated. 

Proof. The corollary follows from the fact that ha is strongly Rayleigh. Such 
measures are negatively associated, which was proved in Theorem 4.9] □ 
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